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ABSTRACT 

We show how one can estimate the multipole moments of the space-time, assuming 
that the quasi-periodic modulations of the X-ray flux (quasi-periodic oscillations), 
observed from accreting neutron stars or black holes, are due to orbital and precession 
frequencies (relativistic precession model). The precession frequencies f2 p and £l z can 
be expressed as expansions on the orbital frequency f2, in which the moments enter the 
coefficients in a prescribed form. Thus, observations can be fitted to these expressions 
in order to evaluate the moments. If the compact object is a neutron star, constrains 
can be imposed on the equation of state. The same analysis can be used for black 
holes as a test for the validity of the no-hair theorem. Alternatively, instead of fitting 
for the moments, observations can be matched to frequencies calculated from analytic 
models that are produced so as to correspond to realistic neutron stars described by 
various equations of state. Observations can thus be used to constrain the equation of 
state and possibly other physical parameters (mass, rotation, quadrupole, etc.) Some 
distinctive features of the frequencies, which become evident by using the analytic 
models, arc discussed. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) open a window to some 
extreme physics. On the one hand the accretion in these 
systems takes place in the region around a compact object 
(such as a black hole or a neutron star) where there are 
strong gravity effects. On the other hand, when the accre- 
tion takes place around a neutron star, the properties of the 
orbital motion depend on the structure of the neutron star 
itself. Thus studying these systems could provide insight on 
the properties of gravity in regions where it is yet untested 
and at the same time help us explore the properties of mat- 
ter in such extreme conditions as the ones in the interior of 
neutron stars. 

One of the properties that LMXBs often display, are 
quasi-periodic oscillations (QPOs) in their X-ray flux. This 
variability appears in a variety of behaviors and is usually 
organized in various groupings, such as the high-frequency 
phenomena, the hectohertz QPOs, or the low-frequency 
complex. The aforementioned categorization is based not 
only on the frequency range but also on the characteristics 
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that the observed frequencies display which depend on the 
type and the state of the source. For example one can ob- 
serve, (i) in the low-frequency complex low-frequency QPOs 
with v values in the range of 0.01 up to 100 Hz or (ii) in the 
case of high-frequency phenomena, neutron star kHz QPOs 
(often observed in pairs), called tw in kHz QPOs , which have 
been observed up to v = 1258 Hz (|jonker et all (|2007h V On 
the other hand there are high-frequency QPOs which are 
categorized as black hole high-frequency QPOs and appear 
to be different than their neutron star counterparts. Finally, 
the QPOs, apart from the previously mentioned differentia- 
tion exhibit also correlations between the various frequency 
components, bot h in neutron stars as well a s in b lack holes 
(for a review see lLambl (120031 ) Ivan der Klisl l|2006h ). 

To explain the QPOs, various types of mechanisms have 
been proposed. These are: (i) the Beat-frequency models, 
where one assumes that there is some "beating" of an or- 
bital frequency by the spin frequency of the central object, 
(ii) the relativistic precession models, where the QPOs are 
associated with the orbital motion and the periastron or 
nodal precession of a particular orbit, (iii) the relativistic 
resonance models, where a type of resonance between the 
orbital and the epicyclic frequencies is assumed wherever 
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they have simple integer ratios, and finally (iv) the preferred 
radii models, where some mechanism chooses a particular 
radius. These models generally assume the geodesic or al- 
most geodesic orbits of the fluid elements in the accretion 
disc to be the sou rce of the observed frequencies (see e.g. 
Ivan der kuI Hooi)), while there are also models in which 
the frequencies are produced from o scillat ory modes of the 
entire disc (see e.g. iRezzolla et alj (|l995l )). In one way or 
another all of these models use the properties of the orbits 
around the compact object onto which the accretion takes 
place. In our discussion we will refer to the models that 
assume that the QPOs are caused by the frequencies associ- 
ated with the orbital motion of the material in the acc retion 
disc s uch as the relativistic precession models (see IStellal 
(2001)). but our analysis has relevance to the other models 
as well. 

In the relativistic precession models it is assumed that 
there is some structure or inhomogeneity in the accretion 
disc (i.e. a hot blob) that follows a keplerian or geodesic orbit 
which is almost circular and almost at the equatorial plane. 
Such an orbit has an orbital frequency Q, and precession 
frequencies, fl p of radial perturbations and f2 z of azimuthal 
(vertical) perturbations. We note here that in the following, 
while we will assume that such a structure is incorporated 
in an accretion disc, the orbit of the structure is such that 
it can be considered to be geodesic, i.e. all forces other than 
gravity can be ignored. This can be considered to be a good 
approximation in the case of radiatively efficient thin discs. 

The most commonly used assumption for describing the 
aforementioned frequencies is that the background space- 
time is either Schwarzschild or Kerr or that it is well ap- 
proximated by one or the other. A more general assumption 
for the background space-time would be to assume a sta- 
tio nary and axi s ymm etric space-time as the one given by 
the IPapapetroul ( 19531 ) line element, 

ds 2 = -f (dt - udcf)f + r 1 [e 27 {dp 2 + dz 2 ) + p 2 dcj> 2 ] . (1) 

Such a space-time ca n be de scribed with the help of 
the Ernst potential (|Ernstl (|l968al lbl)) and the met- 
ric functions can be asso ci ated to the m u ltipole mo- 
ments defined by iGerocbl dl970l ); iHansenl (| 19741 ) and 
iFodor. Honselaers fc Periesl (1 19891 ). Thus one can express 
the rotational and epicyclic frequencies as functions of the 
multipole mome nts of t he sp ace-time. In particular, it was 
demonstrated bv lRvanl (|l995l ) that the quantities Q. p /Q and 
Q. z /Q. can be expressed as series expansions on Q, where the 
coefficients depend on the multipole moments of the un- 
derlying space-time. Although Ryan's work was intended to 
be applied on the analysis of gravitational waves from the 
adiabatic inspiral of a compact object onto a super-massive 
central object (extreme mass ratio inspirals, EMRIs), the 
same principle can be applied in the case of X-ray observa- 
tions from accretion discs (in LMXBs), assuming that the 
orbits in the disc are almost circular and almost equatorial. 

The expressions given by Ryan are asymptotic ex- 
pansions that do not give a very good description for 
the regions close to the innermost stable circular or- 
bit (ISCO). In that case, one could make use of the 
several proposed analytic space-times that give a faith- 
ful description of the geom e try around neutron star s 
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iTeichmiiller. Frob fc Maucheij |201 ll )) and use them to work 
in the strong gravity region. These analytic solutions are 
constructed by using the first few multipole moments of the 
space-time which are associated to the properties of the com- 
pact object. Thus the assumption of such solutions could 
relate the frequencies of the QPOs, to the moments of the 
compact object and give a more accurate description for the 
QPOs than the typically used geometries of Schwarzschild 
and Kerr. In our a nalysis, we will use the two-so l iton an- 
alytic solution (see iManko. Martin fc Ruizl (1 19951 ); IPappaj 
|2009l )) which is constructed from the first four multipole 
moments of the neutron star, i.e., the mass (M), the angu- 
lar momentum (J), the mass quadrupole (M2) and the spin 
octupole (S3). 

The benefit of this approach is that instead of relying on 
gravitational waves (that are at the moment eluding us) to 
get information for the central compact object, one can rely 
on electromagnetic observations which are currently abun- 
dant. Thus, an algorithm that would produce estimates for 
the multipole moments of the space-time in low mass X-ray 
binaries could be used to probe either the equation of state 
for neutron stars or the no-hair theorem for black holes. Fol- 
lowing this point of view, there is currently a programme put 
forward bv lPsaltd (|2008l ) on testing gravity with observa- 
tions in the electromagnetic spectrum. 

Our aim in this work is to demonstrate how one could 
estimate the first few space-time moments by using the se- 
ries expansions of the frequencies (presented in section 2) 
and fitting them to QPOs (in section 4). We also present 
templates for the relation between the various frequencies 
(in section 5), which are constructed with the help of an an- 
alytic solution of the Einstein field equations which has been 
chosen appropriately so as to describe faithfully the geom- 
etry around neutron stars built from realistic EOSs (briefly 
presented in section 3). Using these templates we show how 
one could constrain the EOS for matter of higher than nu- 
clear densities that exists inside neutron stars (section 6). 



2 THE FREQUENCIES AS FUNCTIONS OF 
THE MULTIPOLE MOMENTS 

In this section we will briefly present the expressions that 
describe the orbital motion of a test particle in an axisym- 
metric space-time near the equatorial plane and outline the 
algorithm that associates the precession frequencies with the 
orbital fre quency and t he multipole moments as it was de- 
veloped bv lRvanl (|l995l ). 

A stationary and axisymmetric space-time can gener- 
ally be written in the form of the Papapetrou line element 
given in Eq. |l}, where /, ui and 7 are functions of the 
Weyl-Papapetrou coordinates (p, z). In that case, the met- 
ric components are: 



Stute fc Camenzind](l2002l);lBerti fc Stergioulas! 
IPachon. Rueda fc Sanabria-Gomej ( 20061 ); Pappas] 



2004); 



2009); 



gtt = -/, gt* = fu, g^ = / 1 p 2 - Jlj 2 , 

r-l 27 

g PP = Qzz = f e 

where the radial coordinate is defined as p 2 = g 2 ^ — gttg<\>$- 
The orbit of a particle is usually calculated from the conser- 
vation of energy, the conservation of the z component of the 
angular momentum (where z is the direction of the symme- 
try axis) and the normalization of the four-velocity. Thus, 
the two conservation equations are: 
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-gtt 



dr 
dt 
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E f dt 
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while the normalization of the four-velocity gives: 



(2) 
(3) 



»(£)'+<» (f)(2) 



d0\ 2 /dp A 2 /d2\ 2 



(4) 



If the orbit of the particle is on the equatorial plane, 
then one could define an effective potential from the equa- 
tion 

(dp\ 2 E 2 g<l>4 , + 2EL z g t(f> + L 2 g tt _ 

where _B and L z are the conserved energy and angular mo- 
mentum per unit rest mass, respectively, and V(p) is the 
effective potential. From this equation one sees that the con- 
ditions for circular orbits are, dp/dr = and d 2 p/dr 2 — 0, 
which are equivalent to the conditions for a local minimum 
of the effective potential: V(p) = 0, dV(p)/dp = 0. For the 
circular motion, one can also derive the orbital frequency, 



dt 94>4>,p 



Qtt.pQ 



(6) 



In case the orbit is not exactly circular or not exactly 
on the equatorial plane, one can use expressions (1 21 1 D to 
derive the equation for the perturbed circular motion and 
from these the epicyclic or precession frequencies Q p and Q z . 
The expression for the frequencies is 



series with respect to v = (Mfl) 1 ' 3 (where, in the Newto- 
nian limit, v is the particle's orbital velocity), 



"p \ "* D n 

lT=2^ RnV 



Q 



"z \ ^ r, r. 

1T = 2^ Z " V 



(10) 



(11) 



The coefficients of the series, R n and Z n , depend on the 
Geroch-Hansen-Fodor multipole moments of the space-time. 
For example, the first few non-vanishing coefficients are 



Z a = 2 



Si_ 



Za = 



3M 2 _ Si M 2 
2M3' 6 ~ V M* +6 AP 



where M and M2 are the mass and the quadrupole mo- 
ment, respectively, and Si is the angular momentum. These 
expre s sions extended up to terms of 0(v 1{ ) can be found in 
lRvanl(|l995l ). 

We should point out that in our analysis we are using 
the term "frequency" for the fl values which are actually 
angular velocities while in the literature the observations of 
QPOs are given as "cycles/time" (y values). The two are re- 
lated as fl = 2ttu. Whenever we refer to the v frequencies it 
will be clearly stated so. All quantities are calculated in ge- 
ometric units (km) except for the frequencies which are usu- 
ally given in kHz. The conversion from km -1 to kHz is per- 
formed by multiplying by 299.79, i.e. lkm -1 = 299.79 kHz. 
In these units 1 solar mass (Mq) is 1.477 km. 



{gu +gt^f 



g±±\ 
p 2 ) 



( 



(7) 



Hgt<t> + g^f ( ^ 



-, 1/2 



where if we set a — > p, we obtain the periastron precession 
and if we set a — > z, we get the nodal precession. Thus, we 
have the expressions for the orbital frequency and the preces- 
sion frequencies of an almost circular and almost equatorial 
orbit. 

These expressions depend on the metric functions and 
their first and second derivatives with respect to the coor- 
dinates p and z , evalu ated on the equatorial plane. As de- 
scribed in lRvanl(|l995h . these functions can be expressed as 
power series in 1/p. So, by replacing these expressions in the 
definition of the orbital frequency we obtain 

-l/2s 



fi = (A//p 3 ) 1/2 (l + 



series m p 



This series can be inverted to take the form 
1/p = (Q. 2 /M) 1/S (l + series in 1/3 ). 



(8) 



(9) 



If we use the above expression to eliminate the 1/p de- 
pendence of the precession frequencies, we end up with a 



3 THE SPACE-TIME AROUND NEUTRON 
STARS 

In this section we will present the method that we have 
used to model the space-time around neutron stars. First 
we briefly discuss the procedure used for generating analytic 
space-times and then we talk about the choice of specific 
neutron star models. 



3.1 The analytic space-time 

As we have said previously, a stationary and axially sym- 
metric space-time can be described by the line element ([I}. 
In that case the Einstein field equations in vacuum reduce 
to the Ernst equation (see [Ernst] (|l968al )). 



Re(S)V 2 S = V£ • VS, 



(12) 



where the Ernst potential £, is a complex function of the 
metric functions. 

A general procedure for generatin g solutions of th e 
Ernst equation has been developed by Sibgatulin (1 1991T) ; 
iManko fc Sibgatulinl (Il993f): iRuiz. Manko fc Martini (fl995h 
and IManko. Martin fc Ruid (|l995l ). The solution of the 
Ernst equation is produced from a choice of the Ernst po- 
tential along the axis of symmetry of the form 



S(p = 0,z) = e(z) = 



R(z)> 



(13) 
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where the functions P(z), R(z) axe polynomials of order n 
with complex coefficients in general. In Ruiz et al. (1995), 
one can find a detailed description of the procedure for gen- 
erating the solutions by using the parameters of the polyno- 
mials in the Ernst potential. It is shown in Ruiz et al. (1995) 
that the metric functions are finally expressed in terms of 
some determinants. 

The vacuum two - solito n solution (constructed by 
iManko. Martin fc Ruiz] (|l995l )) is a special case of the pre- 
vious general axisy mmetric solution that is obtained from 
the ansatz (see also ISotiriou fc Pappasl (|2005l )) 



e{z) = 



(z — M — ia)(z + ib) - k 



(14) 



(z + M - ia)(z + ib)-k' 

where all the parameters are real, while the parameters M 
and a are the mass and the reduced angular momentum J/M 
respectively. For this analytic solution, the first five mass and 
mass-current moments of the corresponding space-time are 



M = M, Mi = 

^a 4 — (3a 2 — 

0, J = aM, S 2 = 0, 
-(a 3 - (2a-b)k)M, S 4 



Ah 

So 

S 3 



M 2 = -(a -k)M, 
2ab + b 2 )k + k 2 + - (2a 



0. 



M 3 = 0, 



k)M 2 ^j 



M 



(15) 



As it can be seen, the first four moments are linearly related 
to the parameters of the Ernst potential and thus if they 
are specified, one can generate an analytic space-time with 
exactly these first four moments. Specifically, the parameter 
k is uniquely fixed by the quadrupole M 2 and the parameter 
b is uniquely fixed by the spin octupole 53, while M and a axe 
the mass and the angular momentum per mass, respectively. 

The two-soliton solution has been thoroughly studied 
and has been shown to be a faithful description of the 
exterior of neutron stars and captures all the cha racteris- 
tics of the corresponding numerical space-times fsee lPappasI 
(HH), and for more details ?). 

3.2 Neutron star models 

In order to construct the analytic space-time exterior to a 
compact object, one has to choose the appropriate mul- 
tipole moments. For neutron stars, these moments can 
be read from numerical models constructed by assum- 
ing realistic equations of state. There are several schemes 
developed for numerica l ly int egrating stellar models (see 



IStergioulas fc Friedman (1 19951). and for a n extended list of 
numerical schemes see IStergioulas! (|2003T l). We have used 
Stergioulas's RNS code for the construction of the models. 

The parameters used to construct the analytic space- 
time models, i.e., M, a, k and b, are evaluated from the first 
four multipole moments (M, J, M 2 , S3) as given by equa- 
tions (|15[) . In order to cover more space on the "neutron-star 
parameter space" , these moments are obtained from numeri- 
cally calculated neutron stars with EOSs of varying stiffness. 
For that purpose we have chosen AU as a typical soft EOS, 
FPS as a representative moderate stiff EOS and L to de- 
scribe stiff EOS. After the numerical models are computed, 
the multipole mo ments are evaluated according to the algo- 
rithm described in lPappas fc Apostolatosl ()2012l ) since there 
is a systematic discrepancy in the methods used up to now 
in the literature. 



Table 1. Frequencies calculated for the models constructed with 
the L EOS. The mass and the radii (circumferential radii at the 
equator and at ISCO) are given in km. The frequencies are given 
in kHz. The frequencies are computed at the ISCO except for the 
cases that the ISCO falls bellow the surface of the star; then the 
corresponding frequencies are computed on the surface and are 
indicated by an asterisk. We note that the frequencies shown are 
angular frequencies, i.e. (I = 2irv. 



M 


J 


Rs 


-^isco 




n p 




(km) 




(km) 


(km) 


(kHz) 


(kHz) 


(kHz) 


2.080 


0. 


14.96 


12.48 


7.47* 


4.42* 


0* 


2.071 


0.194 


15.18 


11.76 


7.26* 


3.78* 


0.062* 


2.075 


0.324 


15.61 


11.98 


7.00* 


3.47* 


0.026* 


2.080 


0.417 


16.06 


12.32 


6.74* 


3.25* 


-0.019* 


2.083 


0.483 


16.49 


12.63 


6.51* 


3.06* 


-0.054* 


2.087 


0.543 


16.97 


12.94 


6.26* 


2.87* 


-0.085* 


2.090 


0.598 


17.51 


13.23 


5.99* 


2.66* 


-0.108* 


2.095 


0.650 


18.18 


13.53 


5.68* 


2.42* 


-0.122* 


2.096 


0.690 


18.90 


13.77 


5.37* 


2.18* 


-0.125* 


2.097 


0.698 


19.09 


13.82 


5.29* 


2.12* 


-0.124* 


3.995 


0. 


13.72 


23.97 


5.10 


5.10 





4.012 


0.178 


14.23 


21.83 


5.81 


5.81 


0.144 


4.029 


0.280 


14.69 


20.78 


6.24 


6.24 


0.240 


4.051 


0.375 


15.24 


19.98 


6.62 


6.62 


0.327 


4.074 


0.458 


15.87 


19.45 


6.91 


6.91 


0.390 


4.098 


0.528 


16.53 


19.14 


7.11 


7.11 


0.426 


4.120 


0.588 


17.24 


18.99 


7.24 


7.24 


0.438 


4.139 


0.635 


17.95 


18.95 


7.30 


7.30 


0.434 


4.160 


0.682 


18.93 


18.98 


7.34 


7.34 


0.417 


4.167 


0.700 


19.45 


19.01 


7.35 


7.35 


0.407 


4.321 


0.478 


14.90 


19.77 


6.89 


6.89 


0.496 


4.321 


0.479 


15.01 


19.80 


6.88 


6.88 


0.490 


4.324 


0.483 


15.16 


19.80 


6.88 


6.88 


0.489 


4.325 


0.489 


15.29 


19.79 


6.89 


6.89 


0.490 


4.333 


0.505 


15.55 


19.73 


6.93 


6.93 


0.498 


4.355 


0.555 


16.29 


19.53 


7.07 


7.07 


0.525 


4.377 


0.602 


16.99 


19.40 


7.18 


7.18 


0.541 


4.396 


0.641 


17.68 


19.35 


7.25 


7.25 


0.544 


4.418 


0.684 


18.67 


19.35 


7.30 


7.30 


0.535 


4.420 


0.686 


18.74 


19.36 


7.30 


7.30 


0.534 



For the specifics of the v a rious models chosen, we have 
followed iBerti fc Stergioulasl |2004l ). Thus, we have con- 
structed the same constant rest-mass s equences as the ones 
presented in lBerti fc Stergioula s (2004) for the correspond- 
ing equations of state. That is, every sequence is comprised 
of models with the same baryon number but different ro- 
tation. Specifically for every equation of state, 3 sequences 
of 10 models with varying rotations are constructed, which 
correspond to: 

(i) a sequence terminating to a neutron star of I.4M0 in 
the non-rotating limit, 

(ii) a sequence terminating at the maximum-mass neu- 
tron star in the non-rotating limit, 

(iii) a supermassive sequence that does not terminate at 
a non-rotating model. 

All the sequences end at the mass-shedding limit on the 
side of fast rotation, i.e., at the limit where the angular ve- 
locity of a particle at the equator is equal to the Keplerian 
velocity at that radius. These sequences are the so called 
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Table 2. Fitting results for FPS EOS. The relative differences between the fitted moments and the actual moments of the space-time 
used are given as a percentage. 



Q: 1-3 kHz Q: 0.3-3 kHz Q: 0.1-3 kHz 
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J M J M 2 M J A/ 2 M J A/ 2 
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evolutionary sequences in the sense that an isolated neutron 
star (i.e. with constant number of baryons) that is losing 
angular momentum will evolve along such a sequence. In 
the case that there is accretion onto the neutron star, the 
baryon number is not constant and the exact evolution of the 
neutron star depends on several parameters. For a LMXB 
though, where the mass of the companion is < IMq, the 
overall mass transfer is probably a small fraction of a solar 
mass and thus a neutron star (with mass ^ 1.4Mq) would 
evolve following approximately such a sequence. The choice 
of these particular sequences was made so that the param- 
eter space of the equilibrium neutron star models would 
be fully sampled (for details on neutron star models and 
the ev olutionary sequences, see lCook. Shapiro fc Teukolskvl 
|l99i)). 

Once all the necessary parameters are at hand the var- 
ious frequencies can be easily calculated. As an example, 
we present in Table 1 the radii of the ISCO, the orbital 
frequency SI and the two precession frequencies Q p and Q, z 
at that radius, for the 3 sequences constructed with L EOS. 
We also present the masses, the spin parameters (j = J /Ad 2 ) 
and the surface radii of the corresponding models. In some 
models the surface of the star overcomes the radius of the 
marginally stable orbit and in these cases the respective fre- 



quencies are computed on the radius of the surface. One no- 
tices that the periastron precession frequency Q p coincides 
with the orbital frequency Q at the ISCO and that for the 
non-rotating models, the nodal precession frequency Q z is 
zero, as expected. The ISCO radii are calculated from the 
analytic metric, i.e., the two-soliton metric, and thus deviate 
from the ones given by Stergioulas's code by at most 3 per 
cent jPappaj l|2009l ); ?). 



4 FITTING THE FREQUENCIES USING THE 
ASYMPTOTIC EXPRESSIONS 

In several occasions, there have been observed twin kHz 
QPOs with the upper QPO being in a relativ ely low fre- 
quency range (see eg. iBoutloukos et al.l (|200fj| )) and could 
thus be considered to come from a region of the accretion 
disc further than the ISCO (depending on the mass of the 
central object). In these cases one could directly apply the 
asymptotic expression presented by Ryan and try to fit for 
the multipole moments. 

In order to explore whether this method could provide 
useful measurements of the first moments, we will use simu- 
lated data [pairs of (fi, S7 P )] for the "observed" frequencies 
of QPOs constructed by assuming the two-soliton analytic 
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solution. For every numerical model of those discussed in the 
previous section with multipole moments M, J, M2, S3, a 
corresponding two-soliton space-time is constructed and the 
orbital and precession frequencies for various radii are cal- 
culated. From these frequencies, we select the (f2, Q p ) pairs 
and plot the ratio Q p /Q as a function of Q. We then fit these 
data with the expression (|10[) and calculate the coefficients. 
From these coefficients, the fitted moments are calculated 
and compared to the ones used initially to construct the 
two-soliton space-time. 

For the fits we will use data points which correspond 
to orbital frequencies lower than 3 kHz, which depending 
on the specifics of the model (mass, rotation and higher 
moments) correspond to distances greater than 1.4-2.5 Risco- 
Specifically we will use three frequency ranges for the fits, 
i.e., from 3 kHz down to 1 kHz (which is between 1.4-5 Risco 
depending on the model), 3-0.3 kHz (1.4-10 Risco) and 3-0.1 
kHz (1.4-20 Risco)- As it will be discussed in the following, 
the fit in the moments improves as we increase the frequency 
range, so these ranges have been chosen to demonstrate that 
effect. 

QPOs are observed in the s pectrum of X-ray photons 
in the range of 2-60 keV (see e.g. Ivan der Klisl (|2006l ); lLambl 
(|2003h ). i.e., within a couple of orders of magnitude in en- 
ergy. Generally in an accretion disc, the temperature of the 
disc decreases with the distance as oc {r/Risc o y zli (see 
e.g. iKrolikl (jl999T )l which means that in the range between 
~ (1 — 20) Risco there would be a decrease in emitted pho- 
tons energy of an order of magnitude. Thus, although the 
20 -R^co is not a limiting distance in any sense for the disc or 
the QPOs, it is a reasonable choice for a final fitting range 
for our purpose. 

The results of the fits for the models that are based 
on the FPS EOS are presented in Table 2. These results 
are typical for all three EOSs. One sees from Table 2 that 
in the range of frequencies between 1 and 3 kHz (which 
correspond to distances up to 5Ri aC o), the mass can be very 
accurately constrained. The same applies for the angular 
momentum for most of the models and especially for those 
with faster rotation. The quadrupole is generally not very 
well constrained but for the cases of fast rotation, one could 
have accuracies better than 30 per cent. As the fitting range 
increases, the accuracy of the fits improves and becomes 
quite good for the quadrupole as well. One should also note 
that the accuracy of a fit deteriorates with increasing mass. 
That is because the higher the mass of the model, the closer 
we get to the ISCO for the selected frequency ranges. 



5 FITTING THE FREQUENCIES BY USING 
TEMPLATES 

Since Ryan's expressions are asymptotic expansions in v — 
{MQ) 1 ^ , one should not expect them to be very precise in 
giving the relationship between the frequencies in the region 
of the space-time near the marginally stable orbit, which cor- 
responds to Q. frequencies, greater than about 3-4 kHz. For 
this reason, it would be useful to have templates of frequen- 
cies, for different models constructed with different EOSs, 
on which to project the observations. In this section we will 
present such templates using the three EOSs that we have 
already discussed. 




fl(kHz) 



n (kHz) 



C 0.10 

a 





n (kHz) 



£1 (kHz) 



Figure 1. Plot of frequencies for the three sequences of mod- 
els for the L EOS. Plots (a-d) are for sequence (i), (e-h) are for 
sequence (ii) and (i-1) are for sequence (iii). All curves are ter- 
minated at the ISCO. The black triangles represent the surface 
of the neutron star for the particular model, if the ISCO is un- 
der the surface. The dashed lines correspond to the frequencies 
for Kerr space-times of the same mass and rotation parameter as 
the corresponding neutron star models. The dashed-dotted line 
denotes the maximum QPO frequency measured up to now. 



For every EOS, we produce by means of the two-soliton 
metric, the three frequencies Q, Q p and Q z and the two os- 
cillation frequencies, n p — Q — f2 p for the radial perturbation 
and k z = f2 — Q z for the vertical perturbation. The preces- 
sion frequencies, actually the ratios Sl p /Q, Q z /Q, and oscil- 
lation frequencies are plotted against the orbital frequency 
$1 for the three sequences of neutron star models. In Fig. 
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1 we present the plots for the L EOS. The behavior of the 
frequencies is similar for all the equations of state. 

As it can be seen from the templates, in the case of the 
orbital frequency and the periastron precession frequency, 
as the orbit approaches the marginally stable orbit the two 
frequencies become equal. That is to be expected, since the 
precession frequency is the difference between the orbital fre- 
quency and the frequency of the radial perturbation, which 
goes to zero at the ISCO (the change of sign from positive 
to negative of the square of the frequency of the radial per- 
turbation marks the onset of the radial instability). That 
could be used as a criterion of whether a pair of observed 
kHz QPOs are the orbital and the precession frequencies or 
not. If they are, then one would expect that the evolution of 
the two frequencies should bring them together on the ISCO 
before they disappear. 

Another important feature that arises is that the pre- 
cession frequencies behave differently than the correspond- 
ing Kerr frequencies. More specifically, in some cases the 
nodal precession at high orbital frequencies (which corre- 
spond to the inner part of the disc near the marginally sta- 
ble orbit) becomes zero. That is because the vertical per- 
turbation of the orbit has a frequency equal to the orbital 
one. That never happens in the Kerr geometry (see Fig. 1) 
and the effect is due to the difference between the multi- 
pole moments of a realistic neutron star and those of the 
Kerr space-ti me from the quadrupo l e and higher. It has 
been shown by lLaarakkers &: Poissor] (| 19991 ) that while the 
Kerr quadrupole is equal to minus the square of the spin pa- 
rameter, q = — j 2 , where q = M2/M 3 , for realistic neutron 
stars the quadrupole becomes proportional to that of Kerr's 
with a proportionality constant greater than 1 and increas- 
ing with increasing stiffness of the EOS (it also decreases at 
higher masses). Similar beh avior is observed for the spin o c- 
tupole as well (discussed in lPappas fc Apostolatosl |2012l )'). 
From that perspective a Kerr black hole behaves as an object 
with a super-soft EOS. 

Going further towards the inner regions of the disc 
(closer to the ISCO), the frequency of the vertical pertur- 
bation becomes greater than the orbital frequency up until 
the ISCO. This is evident in the plot of fl z /fl versus SI in 
the case of the L EOS [for the models of the sequence (i) of 
the corresponding neutron stars]. That effect would produce 
a distinctive behavior in the observed QPOs that could be 
used to select an equation of state, since whether or not such 
an effect is present or specifically the frequency range where 
one would expect to see it, depends on the stiffness of the 
EOS of the neutron star. 

We can see similar behavior in the corresponding plots 
for the FPS and AU EOSs, with the difference that the 
vertical oscillation frequency does not always become greater 
than the orbital frequency before the ISCO is reached, as 
can be seen in Figs 2 and 3. The distinctive drop of the 
ratio n z /fi though remains. 

In the plots we also indicate by a small triangle the 
surface of the star whenever it is at a radius that is further 
than the ISCO. Generally the ISCO is swallowed by the star 
for models of the first sequence for every EOS. For these 
models, the decrease of Sl z /Q is also more prominent, thus 
it becomes apparent even if the ISCO is swallowed by the 
star. As it can be seen for the case of L EOS, it might become 
negative even before reaching the star's surface. 
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Figure 2. Cl z /Cl as a function of fi for the FPS EOS. The spin 
parameters j = J/M 2 are also indicated for every model. 
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Figure 3. Same as Fig. 2 for the AU EOS. 

Finally, a dashed-dotted line marks the position of the 
maximum measured QPO frequency (y — 1258 Hz) which 
corresponds to SI ~ 7.9 kHz. The plots show something that 
has been discussed in the literature, namely the fact that this 
particular frequency is incompatible with the L equation of 
state, since there is no model that could produce such a 
high frequency either because of the position of the ISCO or 
because of the position of the surface. 



6 DISTINGUISHING DIFFERENT 
EQUATIONS OF STATE OF NEUTRON STARS 

In this section, we will use the templates presented in the 
previous section and we will attempt to draw some con- 
clusions about the EOS from observations. For illustra- 
tive purposes we will use observations from two sources as 
well, i.e., Scorpius X-l (Sco X-l) represented by the filled 
squares and Circinus X-l (Circ X-l) represented by the 
emp ty squares, as obtained from Ivan der Klis et ail (|l997f ) 
and lBoutloukos et all (|200rj) . 

Once we have the frequency templates for the various 
models that correspond to different EOSs, we can compare 
them to observations of the orbital frequency and of the 
periastron precession frequency, i.e., the twin kHz QPOs. 
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n (kHz) 

Figure 4. QPOs observations compared with the templates for 
the equations of state (a) L, (b) FPS, and (c) AU. The black 
triangles indicate the surface whenever it is present. The filled 
squares are observations related to the source Sco X-l and the 
empty squares are observations of Circ X-l. 



Such a comparison is shown in Fig. 4, where we have plotted 
the curves of Sl p /Q as a function of Q for all the models of 
the three EOSs. The plots also show the observations of Sco 
X-l and Circ X-l. An appropriate selection of observations 
of several sources compared with the frequency templates 
could select or exclude EOSs. For example the two sources 
shown here, appear not to fit very well with the FPS EOS. 
Furthermore, if one had good observations on one source, 
then one could select a specific model and thus estimate the 
mass and the spin of the neutron star, as well as its higher 




500 2000 



V U (Hz) 




V„ (Hz) 

Figure 5. Simulated QPO frequencies related to the periastron 
and nodal precession (y p , v z ) as functions of the orbital frequency 
{v u ) for the models with equations of state (a) L, (b) FPS, and 
(c) AU. Since Q z could be either positive or negative, the u z is 
simply |f2 z /27r|. The characteristic features that appear on the 
nodal precession frequencies of the neutron stars do not appear 
in the corresponding frequencies for a Kerr background. 



moments. That would offer the opportunity for investigating 
the interior structure of a neutron star. 

In Fig. 5 we have plotted for the different EOSs what 
would be the observed frequencies v of the various mod- 
els, assuming that the upper kHz QPOs correspond to the 
orbital frequency, the lower kHz QPOs correspond to the 
periastron precession and the low-frequency QPOs corre- 
spond to the nodal precession. The plots also indicate with 
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a dashed curve the surface for the models that have their 
ISCO inside the star. This type of plots are often used to 
present data from various sources. For illustrative purposes 
we also show with black squares the kHz QPOs of Sco X- 
1 and with empty squares the kHz (upper grouping) and 
low- frequency (lower grouping) QPOs of Circ X-l. Although 
such a plot "masks" the differences of each source when we 
are considering kHz QPOs, it makes evident the behavior of 
the frequency of the nodal precession (low- frequency QPOs), 
that we discussed in the previous section. The turning of the 
low-frequency QPO to a horizontal path and the subsequent 
drop of the curve as the upper kHz QPO increases, that is 
evident in the L EOS plots, would be a clear signature that 
the observed frequencies are indeed related to orbital mo- 
tions. The range of frequencies that such a behavior would 
be observed would also distinguish between different EOSs, 
as the three plots show, while with sufficient resolution one 
could identify a specific model from the frequency of that 
"drop" . 

It should also be pointed out that the three different 
sequences of neutron star models that we have used group 
together in a similar fashion. The first sequence that cor- 
responds to neutron stars that are not near the maximum 
stable mass at the non rotating limit (M — IAMq) exhibits 
lower low-frequency QPOs that present the characteristic 
"drop" discussed previously, while the other two sequences 
(the one at the maximum stable mass in the non-rotating 
limit and the one with no stable non-rotating limit) group 
together to higher low-frequencies and do not show the afore- 
mentioned "drop" . Differences between the three sequences 
for the kHz QPOs are not that prominent, but depending 
on the EOS one could detect a small gap between the first 
sequence and the other two. 



7 CONCLUSIONS 

We have investigated the merits of applying some new tools 
in the study of compact objects from the observation of 
QPOs in the X-ray flux of accreting systems. 

We have shown that Ryan's asymptotic expressions for 
the frequencies, combined with the appropriate observa- 
tions, can be useful in constraining the first three multipole 
moments of the compact object. That would offer us the op- 
portunity to probe the structure of neutron stars, as well as 
the possibility of testing the validity of the no-hair theorem 
for black holes. 

We have also demonstrated the potential benefits of us- 
ing analytic space-times that are constructed so as to be 
faithful to the geometry around realistic neutron star mod- 
els. The comparison of the observed QPOs to the ones cal- 
culated with such space-times, could potentially select or 
exclude EOSs for neutron stars or could even select specific 
neutron star models for specific sources. We have also shown 
that there are characteristic properties of the frequencies of 
the orbits in realistic neutron star space-times that the Kerr 
and Schwarzschild geometries cannot capture. Specifically 
the way that the nodal precession varies as one goes from 
lower to higher orbital frequencies would be a signature of 
the orbital nature for any observed QPOs that exhibit that 
behavior and could be used to distinguish different EOSs 
and probably even particular models of neutron stars. 



Following this line of research it would be worth pursu- 
ing further, with the help of analytic space-times, the inves- 
tigation of orbits that deviate significantly from the equa- 
torial plane as well as the investigation of possible orbital 
resonances that can be observed for these space-times. The 
use of analytic instead of numerical geometries would greatly 
simplify such an exploration. 
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